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ABSTRACT 

We search for evidence of diffuse Lya emission from extended neutral hydrogen surrounding Lya 
emitting galaxies (LAEs) using deep narrow-band images of the Extended Chandra Deep Field South. 
By stacking the profiles of 187 LAEs aX z = 2.06, 241 LAEs at z = 3.10, and 179 LAEs at z = 3.12, 
and carefully performing low-surface brightness photometry, we obtain mean surface brightness maps 
that reach 9.9, 8.7, and 6.2 x 10~^^ ergs cm~^ s~^ arcsec"^ in the emission line. We undertake a 
thorough investigation of systematic uncertainties in our surface brightness measurements, and find 
that our limits are 5-10 times larger than would be expected from Poisson background fluctuations; 
these uncertainties are often underestimated in the literature. At z ~ 3.1, we flnd marginal evidence 
for extended halos with scale lengths of 5-8 kpc, and demonstrate that sub-samples of galaxies with 
low equivalent widths and brighter continuum magnitudes are more likely to possess such halos. At 
z ~ 2.1, we find no evidence of extended Lya emission down to our detection limits. We compare these 
findings to other measurements in the literature, and discuss possible instrumental and astrophysical 
reasons for the discrepancies, including possible evolution of Lya halo properties from z ^ 3.1 to 
z - 2.1. 

Subject headings: galaxies: evolution - galaxies: high-redshift — galaxies: halos - galaxies: structure 
- cosmology: observations 



L INTRODUCTION 

The star formation rates of high and intermediate red- 
shift galaxies depend on the complex interplay between 
accretion and outflow. While infalling gas f rom the cos- 
mic w eb provides fuel for starbursts (e.g., iDekel et al.l 
|2009( ). the winds of young stars and supernovae cause 
feedback and dampe n the intensity of the accretion (see 
iVeilleux et al.l I2005D . The star formation history of 
a galaxy and its surrounding circum-galactic medium 
(CGM) are thus inextricably intertwined, and any model 
which seeks to reproduce the observed properties of these 
systems must take into account both of these processes 
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(e.g., IDave et al.llMTallbh . 

Observationally, there is strong evidence for the pres- 
ence of galactic-scale outflows around luminous star- 
forming galaxies at high redshift. At z ^ 3, the ve- 
locity offsets between resonance and non-resonance ab- 
sorption features of Lyman- break galaxies (LBGs) (e.g., 
Steidel et aDll996t lAdelberger et al1l2003l : iShaplev et alJ 
20031: ISteidel et al.ll20Tfli r and the relative strengths and 
velocities of multiple-peaked profiles of Lya in emis- 
sion (e.g., Verhamme et al. 2006; Tapken et al. 200^ 
Verhamme et al. 2008; Laursen et al.ll2009t iBarnes et ahl 
2011: Kulas et al. 2012) all point to outflows being com- 
mon in these high luminosity systems. Outflows are 
also comm on in lower- l umin osity objects: the spectral 
analysis of iBerrv et al.l (|2012l ) clearly demonstrates that 
Lya emitting galaxies at z ~ 3 have strong galactic 
winds. However, while the evidence for outflows is com- 
pelling, there is as of yet no direct evidence for the cold 
gas infall that is pr edicted by galaxy formation models 
(ISteidel et aPIMoh . 

Detecting the circum-galactic material around high- 
redshi ft galaxies is a difficult proposition. ISteidel et al.l 
(I20T0I) was able to infer the existence of this component 
out to ~ 120 kpc by using angular pairs of unassociated 
z ~ 2 — 3 star-forming galaxies and looking for Lya and 
Ly/3 in absorption. However, a more direct approach is to 
search for circum-galactic material in emission. Resonant 
scattering can cause the Lya photons produced inside a 
galaxy to scatter hundreds or even thousands of times 
before escaping into intergalactic space. Consequently, 
if a galaxy is surrounded by neutral hydrogen gas, that 
materi al is likely to be illti minated by Lya. In fact, mod- 
els bv iZheng et al.l ()2011| ) suggest that at z = 5.7, dif- 
fuse Lya emission in the CGM may be observable several 
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hundred kpc from the site of its creation. In addition to 
inflows and outflows, a clumpy and inhomogeneous ISM 
may also create halo-like stru ctures. 

The recent observations by iSteidel et al.l (|201lL here- 
after Sll) appear to confirm this prediction. By co- 
adding Keck, Hale, and Subaru narrow-band images of 
92 LBGs between 2.2 < 2 < 3.2, Sll claimed the de- 
tection of scattered Lya from a diffuse halo that ex- 
tends ~ 80 kpc from their mean composite galaxy. The 
azimuthally-averaged surface brightness profile of this 
halo is exponential, with a scale length that is substan- 
tially larger than that of the galaxy's coi itinuum, but 
similar to that measured for Lya lalobs (jSteidel et al.l 
I200a iMatsuda et al.|[200l . This consistency led Sll to 
argue that Lya halos are a generic feature of all high- 
redshift star-forming galaxies, and that these halos have 
self-similar exponential profiles. 

However, the galaxies targeted by Sll were all selected 
via the Lyman-break technique, and are thus amongst 
the brightest and most massive galaxies in the high red- 
shift universe. With stellar masses of ~ 2 x 10^°Mq, 
star formation rates (SFRs) of ^ 30M(t, yr~^, and 
internal extinctio ns o f Ay ^ 1 (jShaplev et al.l 120011 : 
lAdelberger et ar.ll2005D . these objects are in the extreme. 
Most high-z-objects have much lower masses and star 
formation rates, and are not detectable via the Lyman- 
break technique. To test the ubiquity of Lya halos, more 
representative samples of objects are needed. 

Lya emitting galaxies (LAEs) may be a better popula- 
tion for this purpose. Although LAEs can have Lya lu- 
minosities similar to that of Lyman-break galaxies, most 
have star formation rates (SFR ~ 2Mo yr"-'^) and stel- 
lar masses (M < 10^ Mq) that are an order of magni- 
tude less. These are the objects that will likely evolve 
into today L* galaxies like the Milky Way ( [Pirzkal et al.l 
I2OO7I: IGawiser et all [2007t IGuaita et al.Tl2010D . More^ 
over, LAEs are relatively metal-po or and dust-free, with 
internal extinctions Ay < 0.5 (IGawiser et al.l 120071 : 
IGuaita et ahl 120101: lAcauaviva et al.ll201lD . This lack of 
dust makes it more likely that the Lya photons generated 
in the systems' star-forming regions will escape, although 
strong outflows in LBGs may also produce similar escape 
fractions. Finally, spectroscopy has shown that the con- 
tamination rate of photometrically-selec ted LAEs is ex- 
tremely small, less than 2% at z ~ 3.1 (jGronwall et al.l 
I2007t iCiardullo et"aLll2012D . Thus, large samples of Lya 
emitters are relatively strai g htforw ard to obtain. 

Recently, IMatsuda et"all H2012I . hereafter M12) used 
the co-added the images of ~ 2100 LAEs to conclude 
that the properties of Lya halos depend on environment, 
with the emission surrounding LAEs in low density en- 
vironments being much fainter than that for correspond- 
ing objects in overdense regions. However, neither Sll 
or Ml 2 explored the systematic effects which limit the 
depth to which image stacking can be trusted. Invari- 
ably, the true limit of any surface brightness analysis is 
significantly brighter than that which is inferred from 
simple counting statistics. 

The Sll study focused on LBGs at 2.3 < z < 3.1 
while M12 focused on LAEs at z ^ 3.1 in a range of 
galaxy environments. In this paper, we investigate the 
extended Lya halos of field LAEs at similar redshifts. In 
§2, we define our samples using the z ^ 3.1 observations 



levels of « 8 X 10 ergs cm ^ s ^ arcsec ^ at z 
and w 1 X 10~^* ergs cm~^ arcsec"^ at z 



of lGronwall elall (l2007l) andlCiardullo et all (I20TI) , and 
the z ~ 2.1 data of IGuaita et al.l ( 20100 . In §3, we de- 
scribe our stacking technique and present the initial anal- 
ysis of the Lya profiles of our LAE sub-samples. Then 
in §4, we perform a detailed analysis of the various sys- 
tematics which limit the precision of stacking analyses, 
and show that the systematic uncertainties in this pro- 
cedure have generally been underestimated. In §5, we 
present a more detailed stacking analysis, and demon- 
strate that there is mixed evidence for extended Lya- 
emission in our samples, down to a surface brightness 
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show that stacks containing sources that that have rel- 
atively low Lya equivalent widths and/or bright contin- 
uum magnitudes are more likely to show extended Lya 
emission than similar stacks with low continuum lumi- 
nosities or high equivalent widths. In §6, we compare our 
results to previous studies via a series of Monte Carlo 
simulations, and demonstrate that halos such as those 
found by Sll would have easily been detected in our 
survey. In §7, we discuss possible explanations for our 
results, including the extended telescope point-spread- 
function, the limitations of sky subtraction, the intrinsic 
differences between galaxy populations, the effects of en- 
vironment and orientation, and the possible evolution of 
Lya halos between z 3 and z ~ 2. Finally, in §8, we 
summarize our results and discuss the possible ways of 
improving the constraints on Lya halos with upcoming 
surveys and new Lya radiative transfer models. 

For this paper, we assume a ACDM cosmology with 
nx = 0.70, = 0.3 and Hq = 70 km s^^ Mpc^^. With 
these values, 1" — 8.32 kpc at z = 2.1 and 7.63 kpc at 
z = 3.1. 

2. THE SAMPLE OF LYa EMITTERS 

In order to test for the existence of diffuse Lya 
around Lya emitting galaxies, we consider three sta- 
tistically complete sets of LAEs. Each set is located 
in the Extended C handra Deep Field South (ECDF-S; 
iLehmer et al.ll2005l) , each was identified via narrow-band 
imaging with the Mosaic II CCD camera (jMuller et al.l 
IT998) on the CTIO Blanco 4-m telescope, and each de- 
fines an LAE as an emission-line galaxy with a rest-frame 
Lya equivalent width greater than 20 A.^ 

Our first set of LAEs is that found by iGronwall et al.l 
(|2007D using a 50 A full-width-half-maximum (FWHM) 
filter centered 4990 A. This dataset, hereafter called C- 
03, covers the redshift range 3.08 < z < 3.12 and is com- 
plete to a monochromatic Lya emission-line luminosity 
limit of log i = 42.1 ergs s^^. A total of 254 objects were 
identified in this survey, though only 156 are members of 
the statistically complete sample. To date, 67 of these 
LAEs have been spectroscopically confirmed as z ~ 3.1 
star-formi ng galaxies (the firs t 61 confirmations are re- 
ported in IGawiser et al.l l2007| ) . The median seeing for 
this survey was I'.'O. 

The second set of LAEs is that found by 
ICiardullo et all (poTl . using a 57 A FWHM fiffer 
centered at 5010 A. This LAE catalog, hereafter called 
K-03, is complete to logL > 42.3 ergs s~^, and ha s 
considerable overlap with that of lGronwall et all ()2007D . 
However, because the K-03 central wavelength is slightly 
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redder than that of C-03 , the dataset also contains 
68 new objects in the redshift range 3.10 < z < 3.15. 
This K-03 survey discovered 199 LAE candidates; 130 
of these are in the statisticahy complete sample. The 
median seeing for this survey was I'.'l. 

Our third sample of LAEs examines a different red- 
shift regime. By using a 50 A FWHM filter centered 
at 3727 A, iGuaita et all ()2010D were able to identify 
250 L AE candidates i n the redshift range 2.04 < z < 
2.09. ICiardullo eFall (12012) further refined this sam- 
ple, by improving the photometry, removing one proba- 
ble AGN, and using Hubble Space Telescope morpholog- 
ical measurements to exclude 34 likely interlopers. This 
left a final sample consisting of 201 objects, of which 
73 are brighter than a completeness limit of \ogL = 
41.8 ergs s~^. We define this as the 02 LAE sample. 
The median seeing for this survey was 1'.'4. 

Since all three surveys targeted a region of space 
with deep Chandra data, most quasars and other ac- 
tive galactic nuclei have been elimi nated from the sam- 
ples based on thei r X-ray ernission (jLehmer et al.l [20051 : 
iVirani et all 120061: iLuo et al.ll2008D . Similarly, because 
the ECDF-S has been subjected to deep near-UV (1700- 
2800A) surveys by the GALEX and Swift sa tellites 
(jSchiminovich et al.l [2003t iHoverste n et al.ll2009[ ). most 
foreground interlopers have already been identified. (For 
the redshifts under consideration, the rest-frame near- 
UV lies blueward of the Lyman break.) Thus, our sam- 
ples of LAEs should be relatively clean. Of the 73 
2 3.1 LAEs with spectroscopic confirmations, all have 
redshifts that place the galaxy between 3.08 < z < 3.12. 
Similarly, of the 13 z = 2 .1 objects with spectroscopic 
redshifts (jBerrv et al.l20l3) , there are no foreground con- 
taminants. 

Finally, we should note that almost 300 LAEs in the 
E-CDFS (171 at z - 3.1 and 108 at z ~ 2.1) have been 
imaged in the rest-frame ultraviolet with Hubble Space 
Telescope's Adva nced Camer a for S urveys. An analysis 
of these data by iBond et al.l ()2012[ ) has shown that in 
Vgoe , LAEs in the redshift range 2 < z < 3 are extremely 
compact, with median half-light radii less than 0'.'2. In 
fact, at z = 3.1, the largest half-light radius recorded is 
0'.'29 and the largest half-light radius at z = 2.1 is 0'.'43. 
In both cases, this is much smaller than the seeing of the 
ground-based images under consideration (I'.'O and 1'.'4), 
respectively. Consequently, if the Lya emission emanat- 
ing from these objects has the same spatial distribution 
as the UV light, the profiles of these objects should be 
consistent with those of the frames' point spread func- 
tions. 

3. INITIAL STACKING OF LYa EMITTERS 

To search for ext e nded L y a, we began w it h the frames 
of iGronwall et all (I2007D . IGuaita et all (|2010[ ). and 
ICiardullo et al.l (120121)" all of which were processed us- 
ing th e Mosaic II reduction procedures of iG awiser et al.l 
(poM ). First, a small 27" x 27" (101 x 101 pixel) re- 
gion was extracted about the center of each LAE. To 
avoid edge effects, any source close to the CCD bound- 
aries was excluded at this step. This left us with 146, 
124, and 69 sources in the statistically complete C-03, 
K-03, and 02 samples, and 241, 179, and 187 objects 
in the total samples of these surveys. Next, to examine 



the behavior of Lya halos in the more massive LBGs, we 
used the broadband magnitudes of the MUSYC survey 
(|Gawiser et al.|[2006[ ) to define sub-samples of LAEs more 
luminous than Rab = 25.5. This continuum magnitude 
limit and equivalent width properties of this sub-sample 
are similar to the LBGs studied by Sll, although these 
sources were originally selected on the basis of equivalent 
widths, not continuum colors. These sub-samples, which 
we will refer to as "UV bright", contain 28, 19, and 27 
sources in the C-03, K-03, and 02 samples, respectively. 

Finally, to examine whether equivalent width has any 
affect on the properties of Lya halos, we ranked our LAEs 
by their (photometrically-derived) rest-frame equivalent 
width, and created sub-samples of objects in the high- 
est and lowest quartiles of the distribution. We refer to 
these sub-samples as "high EW" and "low EW" sub- 
samples, and each contains 60, 44, and 50 objects in 
the C-03, K-03 and 02 datasets, respectively. For the 
low EW sub-samples, the rest-frame equivalent widths 
range from 20 A to 60 A; for the high-EW sub-samples, 
EWo > 74 A, with most of the objects having rest-frame 
equivalent widths greater than 100 A. We note that since 
the distribution of LAE equivalent width s has (at best) a 
weak dependence on galaxy luminosity ([Gronwall et al.l 
I2007t ICiardullo et"alll2012f ). samples chosen by equiva- 
lent width and UV-brightness are somewhat orthogonal. 
For example, in the C-03 data, only 25% of the low- 
equivalent width LAEs (16 out of 60) are members of the 
UV-bright sample; for the K-03 sample, this fraction is 
~ 35% (16 out of 44), and for 02, the fraction is ~ 44% 
(22 out of 50). A summary of the different sub-samples, 
and other relevant information, is given in Table [T] 

After defining the LAE sub-samples, each object was 
scaled to a common narrow-band flux level and median 
combined, using the imcombine task within IRAEQ, and 
adopting a (3 a) sigma-clipping algorithm with a scheme 
which weighted each source by its total flux. To avoid any 
effects of interpolation, we used only integer pixel shifts 
for each stack; this procedure has the effect of slightly in- 
creasing the radial surface brightness profile. The stack- 
ing procedure was used both on the original narrow-band 
images themselves, and on continuum-subtracted images, 
where the flux level at 5000 A was defi ned using the 
B and V frames of the MUSYC survey (jGawiser et al.l 
I2006f ). For the 02 data, the continuum flux at 3727 A 
was derived from publicly available U and B frames taken 
with the Wide Field Imager (WFI) of the ESO 2.2 m 
telescope (Hildcbrandt et al. 2006) that were reprojected 
to match the MUSYC BVR image of the same field. 
These continuum subtracted images have the advantage 
of decreasing the contamination from nearby continuum 
sources, but come with the penalty of decreased signal- 
to-noise and possible mis-matches in the frames' point- 
spread-functions (PSFs). 

Finally, as a control on our experiments, we also 
analyzed a catalog o f point sources initially found by 
lAltmann et"all ()2006[ ) on the MUSYC broad-band im- 
ages. This catalog, which consists of objects with the 
colors of Galactic K- and M-type stars, was further re- 

^•^ IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement virith 
the National Science Foundation. 
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TABLE 1 

Sub-Sample Properties 



TTlTll o 


TVnmVinr of T , A P.t; 


\^/ 


lV±s_.LJ.l<Xli O^t^lllti 


1 M 1 X A .III vc l^di±^ w i_faiLi\-L 

Exposure Time (hrs) 


EW^ 

T3 o n rrc ( A \ 

rLtinge 


EW^ 
iVieQian j 


C-03 Total 


241 


3.10 


I'.'O 


24.00 


20 - 366 


54 


C-03 Statistically Complete 


146 


3.10 


I'.'O 


24.00 


20 - 366 


51 


C-03 UV Bright 


28 


3.10 


I'.'O 


24.00 


20 - 158 


30 


C-03 High EW 


60 


3.10 


I'.'O 


24.00 


82 - 366 


103 


C-03 Low EW 


60 


3.10 


1"0 


24.00 


20 - 36 


27 


K-03 Total 


179 


3.12 


1"1 


15.67 


20 - 628 


63 


K-03 Statistically Complete 


124 


3.12 


1"1 


15.67 


20 - 628 


60 


K-03 UV Bright 


19 


3.12 


1"1 


15.67 


22 - 49 


23 


K-03 High EW 


44 


3.12 


I'.'l 


15.67 


94 - 628 


125 


K-03 Low EW 


44 


3.12 


I'.'l 


15.67 


20 - 40 


29 


02 Total 


187 


2.07 


1'.'4 


35.75 


20 - 933 


43 


02 Statistically Complete 


69 


2.07 


1'.'4 


35.75 


27 - 933 


70 


02 UV Bright 


27 


2.07 


1'.'4 


35.75 


20 - 60 


26 


02 High EW 


46 


2.07 


1'.'4 


35.75 


74 - 933 


98 


02 Low EW 


46 


2.07 


1'.'4 


35.75 


20 - 28 


23 



^ Equivalent widths arc all rcst-framc, and were determined from photometry. 
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Fig. 1. — On the left is the median-combined statistically com- 
plete sub-sample of 146 C-03 Lya emitters; on the right is the 
median-combined sample of 55 stellar sources from the same image. 
The greyscale stretch is linear for each stack, and ranges from -la 
to 0.15p, where a is the standard deviation of the pixels near the 
corners of each image and p is the peak pixel value in ADU. There 
is no apparent difference between the two stacks. Formally, the 
FWHM of the Lyo emitters is 8% larger than the stellar sources, 
but the uncertainties are of comparable magnitude (5%). 

fined bv lBond et all ([20lTI) using H ST-WFC3 Ea rly Re- 
lease Science images (ERS; IWindhorsielaLllloiil) of the 
GOODS-S region. These objects were re-identified on 
our narrow-band frames, and stacked in a manner iden- 
tical to that for the LAEs. Fifty- five stars formed the 
point-source stack on our C-03 and K-03 images; for 
the 02 data, our control image was created using the 
median of 19 stars. 

Figure [T] compares the median-stacked image of the 
146 z = 3.1 LAEs in the statistically complete C-03 
sub-sample to that of the field's point sources. A careful 
examination of the figure reveals no compelling evidence 
for extended Lya emission. The results for the other 
LAE stacks are similar: there is no clear signature of 
extended Lya emission in any of the LAE sub-samples. 
Table [T] quantifies this result by listing the FWHM of 
a two-dimensional Gaussian fit to each image stack, as 
measured with the f itpsf task within IRAF. The un- 
certainties on the fits for each stack are derived from a 
bootstrap analysis, in which we randomly resampled the 
sources (with replacement) and repeated the imcombine 
and f itpsf steps 30 times for each stack. For com- 
parison, the FWHM of the point-source stacks are also 
shown, and the LAE measurements are tabulated with 
respect to these values. As the table demonstrates, there 
is a suggestion that some of the LAE stacks have broader 



radial profiles than that of the PSF stars. However, the 
uncertainties in this type of analysis are large and no 
obvious trends are apparent between the different sub- 
samples. Moreover, this simple procedure does not fully 
account for the effects of centroiding errors: because the 
LAEs are much fainter than the comparison stars, the 
FWHM of their stacks are expected to be wider due to 
positional uncertainties. A more sophisticated analysis 
is therefore needed to find evidence of Lya halos. 

4. SURFACE BRIGHTNESS LIMITATIONS 

Before proceeding further, we need to examine the re- 
alistic limits of our surface brightness measurements, and 
properly characterize any systematic uncertainties. This 
is not a straightforward process. Surface photometry at 
faint fiux levels has additional uncertainties above and 
beyond those produced by Poisson statistics. Large-scale 
flat-fielding errors, the structure of the point-spread- 
function (PSF) at large radius, and errors in the esti- 
mation of the sky background all introduce systematic 
uncertainties into surface brightness measurements, and 
these errors dominate the error budget. Such uncertain- 
ties are well known in the astronomical literature (e.g., 
Morrison et al." 1994; Zheng et al. 1999; ' Gonzalez et alJ 
2005; Zibctti et al. 2005; Bernstein 2007), but are not 
always applied in practice. 

The dominant systematic uncertainty for ultradeep 
surfac e photometry is usually the large-scale f iat- fie lding 
error (jMorrison et all 119941: iFeldmeier et al.l I2004D . In 
order to empirically determine this error on the C-03, 
K-03, and 2 survey data, w e follo wed t he procedures 
described in IFeldmeier et al.l ()2004[ ) and iRudick et al. 
^2010). In brief, we began by using DAOPHOT (jStetson 
[19871) to find and measure the magnitude of every point 
source on the frame. We then determined a more precise 
small-scale PSF for each image (out to the larger radii of 
~ 5"), while masking out all of the point sources down 
to an estimated level of 0.1/iexp ADU, where texp is the 
exposure time and ADU is the number of counts above 
the sky value. We next identified a saturated star in 
each image, and, after masking all other nearby sources, 
used it to construct an azimuthally averaged PSF out 
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TABLE 2 

Results from Two-Dimensional Gaussian fitting 





C-03 


K-03 


02 


Sample 


FWHM Ratio'' 


FWHM Ratio'' 


FWHM Ratio'' 




(arcsec) 


(arcsec) 


(arcsec) 



Stacked Stars 1.23 ± 0.03 1 1.23 ± 0.04 1 1.40 ± 0.09 1 

Total Sample 1.39 ± 0.04 1.13 ± 0.04 1.28 ± 0.03 1.04 ± 0.04 1.53 ± 0.07 1.09 ± 0.09 

Statistically Complete 1.33 ± 0.05 1.08 ± 0.05 1.32 ± 0.02 1.07 ± 0.03 1.58 ± 0.07 1.13 ± 0.09 

UV Bright 1.36 ± 0.07 1.11 ± 0.06 1.25 ± 0.05 1.01 ± 0.05 1.54 ± 0.17 1.10 ± 0.14 

High Equivalent Width 1.19 ± 0.07 0.97 ± 0.06 1.19 ± 0.04 0.96 ± 0.04 1.48 ± 0.07 1.06 ± 0.08 

Low Equivalent Width 1.35 ± 0.09 1.10 ± 0.08 1.28 ± 0.05 1.04 ± 0.04 1.31 ± 0.08 0.94 ± 0.08 

Difference Total 1.28 ± 0.04 1.04 ± 0.04 1.21 ± 0.23 0.99 ± 0.19 1.09 ± 0.11 0.78 ± 0.09 

Difference Statistically Complete 1.45 ± 0.08 1.18 ± 0.07 1.45 ± 0.18 1.18 ± 0.15 1.62 ± 0.15 1.16 ± 0.13 

' Defined as FWHM (sub-sample) / FWHM (stacked stars), with the uncertainties from each sample added in quadrature. 
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Fig. 2. — The azimuthally averaged surface brightness profile of 
a star constructed for each of our narrow-band images. The in- 
strumental surface brightness scale, fiinst is defined as 1 ADU s"! 
= 25 mag arcsec"^. The red triangles denote the small-scale PSF, 
as measured by bright stars and the DAOPHOT software package. 
The filled squares represent the merger of the DAOPHOT PSF (at 
small radii) with the azimuthally averaged profile of a saturated 
star (at large radii), where the large-scale PSF was normalized to 
the small-scale PSF. For reference, the dashed vertical lines shows 
the radial extent studied by Sll (12") and M12 (8'.'5) in their anal- 
yses. 

to a much larger radius (50" to 60", depending on the 
image and the brig htness of th e chosen star). As was 
found orig inahy by lKm3 (flOTl . the PSF of a ground- 
based telescope extends out to a very large radius, and 
has a dramatic change in slope in the far field. As dis- 
played in Figure [5] once outside of ~ 4" (i.e., ^ 30 kpc 
at the redshift of our objects), the PSF contains flux 
far above that which would be predicted from a simple 
Gaussian extrapolation. Note the difference between the 
PSF of the 02 image and that of the two images near 
5000 A. This is due to the differing patterns of inter- 
nal reflections between the CCD dewar window and the 
narrow-band filter [Slater et al.l l2009). Measurements of 
the large-s cale PSF by oth er authors (see the detailed dis- 
cussion bv lBernsteinH2007l) indicate that the shape of this 
function depends not only on the telescope, instrument, 
and wavelength, but on the observing season as well, as 
observations taken several months apart may have very 



Fig. 3. — The histogram of sky values for each of our block- 
averaged images, binned into intervals of 5.0 xlO"" ADU/s bin~^. 
Ideally, the sky values should all equal zero. However, due to large- 
scale flat-fielding errors, and the faint unmasked wings of stars and 
galaxies, there is usually a dispersion about zero in observed data. 
This distribution determines the limiting surface brightness of our 
analysis. 

different power-law-slopes at large radii. Consequently, 
any analysis which relies on subtracting the large-scale 
PSFs of different telescopes and/or instruments must be 
treated with suspicion. To be conservative, the remain- 
der of this paper will focus solely our our narrow-band 
data alone; no attempt will be made to interpret our 
continuum-subtracted images. 

With the large-scale PSF determined, we re-masked 
all of the point sources in the frame, again down to a 
level of 0.1/fcxp ADU. We then masked out all of the 
frame's non-poi nt sources, using the SEx tractor software 
package y2.2.2 (iBertin fc Arnoutsl [ l996[) and the proce- 
dures of lFeldmeier et alj (|2004f ) . Finallv. after applying 
this mask, we block-averaged the remaining pixels into 
49 X 49 pixel (13"1 x 13"1) superbins, and calculated 
the robust mean of ea ch superbin containin g more than 
100 unmasked pixels ([Morrison "eFallliggl . Given the 
large size of these superbins, each of their mean values 
should be equal to the sky subtracted mean of the im- 
age. Instead, as the histograms of Figure [3] illustrate, 
there is a significant amount of dispersion about the sky 
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subtracted mean. This is the effect of large-scale flat- 
fielding errors, and the wings of stars and galaxies that 
remain unmasked despite our procedures. 

We conservatively define the large-scale flat-fielding er- 
ror as the full-width at half-maximum of these distribu- 
tions; these values are given in column 2 of Table [1] As a 
cross-check of our analysis, we also determined the stan- 
dard deviation of the sky background at the edges of 
our LAE and point-source image stacks. These values 
will not be identical to those derived from the superbins 
for two reasons: 1) the effect of scaling each individual 
image to a common flux level, and 2) the differing spa- 
tial distribution of the LAEs compared to the sky his- 
tograms. Nevertheless, the values obtained from the sky 
background of the stacked images are within 50% of the 
adopted values of Table [T] 

The third and fourth columns of Table [T] compare 
the large-scale flat-flclding errors to the median sky val- 
ues of the parent image. These uncertainties are com- 
parable to those found for other de ep irn aging surveys 
(jFeldmeier eFaL|[200l iRudick et al.|[20Tol ). but are not 
as good as the those associated with the original MUSYC 
broad-band survey (see Table 5 of Gawiscr ct al. 2006). 
We attribute this difference to the lower flux values asso- 
ciated with our narrow-band data: although the MUSYC 
images were taken with the same telescope and reduced 
using similar techniques, the increased count rate of the 
background sky depressed the relative importance of ad- 
ditive sky terms. More importantly, our large-scale flat- 
flelding errors are flve to ten times larger that those which 
would be inferred using simple Poisson statistics (see the 
5th column of Table [T]) . This noise floor cannot be re- 
duced by simply including more sources in the stack, or 
binning the overall results: it is intrinsic to the data. The 
final column in Table [T] gives this fundamental limit in 
units of monochromatic surface brightness. 

5. DETAILED STACKING 

With the systematic uncertainties of our datasets bet- 
ter determined, we now begin a more detailed analysis of 
the surface brightness profiles of LAEs. To do this, we 
azimuthally averaged each stacked narrow-band image, 
using a harsh (2.5cr) sigma-clipping algorithm, thereby 
creating a series of radial bins extending 18'.'7 (70 pixels) 
from the source center. Uncertainties for each radial bin 
were calculated using a full error model, including read- 
out noise, photon noise, small-scale flat-flelding error, 
and large-scale flat-fleld ing error, and were de termined 
in a similar manner as in lFeldmeier et ahl (|2Q04 see their 
Appendix for a full description). Next, we determined a 
new local sky value for each stack, by computing the me- 
dian pixel value in a circular annulus between 13'.'4 and 
18'.'7 (50-70 pixels) from the LAE centroid. At z - 2.1, 
this annulus corresponds to a linear distance of between 
111 and 156 kpc; at z ~ 3.1, the range is between 102 and 
142 kpc. In all cases, the local sky was within ~ 1% of 
the original sky value, and so the additional uncertainty 
introduced at this step was small (on average a factor 
of two less than the large-scale flat-fielding error). This 
constant sky value was then subtracted from the radial 
proflles. 

The resulting sky-subtracted radial proflles are plotted 
in Figures |4] through |6l Also shown in these figures are 
the frames' total PSFs (core and large-scale; calculated 



in §4), the radial profiles derived from the stacked set 
of point sources, and the surface brightness limits de- 
fined by the frames' large-scale flat-fielding errors. In 
each flgure, the comparison proflles have been scaled in 
flux to match the three innermost points (0"-0'.'6) of 
the corresponding LAE. We note that the lines repre- 
senting the surface brightness limits are approximate, 
since they represent the mean flat-flelding err or of the 
fleld. Since LAEs exhibit angular clustering (jFranckd 
|2009( ). the large-scale flat-flelding error appropriate to 
the stacked images may be slightly different. Neverthe- 
less, they do illustrate the level at which stacks made 
from the wide- fleld Mosaic imager become suspect. 

To better summarize the profiles displayed in Figures|4] 
through [6l we subtracted the point source profile from 
each LAE sub-sample stack, again using the three inner- 
most points for scaling. Next, following Sll and M12, 
we fit the exponential proflle: 

5(r) = C„exp(-r/r„) (1) 

where C„ is the normalization factor, and r„ is the 
exponential scale length, to the residuals, using a 
Levenberg-Marqua rdt non-linear least squares algorithm 
(jPress et al.l I1992I ). which took into account negative 
fluxes for sky-subtracted bins. The results of these flts 
are given in Table [T] The uncertainties given in the table 
are the internal errors from the fit themselves, and are 
likely to be an underestimate to the true uncertainties. 
For practical reasons, our flt was restricted to the spatial 
scales between 0'.'75 (where seeing effects dominate) to 
S'.'O, where large-scale flatflelding errors begin to make 
all measurements consistent with zero. Moving the inner 
radius outward to match the I'.'O value used by Sll does 
not qualitatively change the results. 

An inspection of Figures |4] through [6] and Table [T] re- 
veals several interesting patterns. At 2 3, there is some 
evidence for the existence of extended Lya halos in some 
of the LAE sub-samples. Speciflcally, there is marginal 
evidence for UV continuum bright LAEs to be slightly 
more extended than the general LAE population, and 
for low-equivalent width objects to have larger halos than 
their high-equivalent width counterparts. Since the effect 
is present on stacks made from both the C-03 and K-03 
frames, this suggests that the effect is real, and that there 
is a genuine difference between LAEs with strong and 
weak Lya emission. The sub-samples that show possible 
evidence for extended Lya halos are plotted in Figure [7] 
in residual form, along with the best-fltting exponential 
proflle. For comparison purposes, the "LAE only" expo- 
nential fit of Sll and the "—1 < 51^^ < 0.5 exponential 
flt of Ml 2 are overplotted. 

On the other hand, z 2 LAEs show no evidence 
for this dichotomy, as all the proflles displayed in Fig- 
ure [6] are all consistent with point-source emission. At 
face value, this is peculiar, since the effects of (1 -I- z)'^ 
cosmological surface brightness dimming should cause 
2; ~ 3 halos to be more than ~ 3 times fainter than 
their z ^ 2 counterparts. While the poorer seeing and 
slightly brighter sky of the 02 stacks may partly compen- 
sate for this effect, the lack of extended emission around 
the z 2 objects may also indicate differences in the 
populations of LAEs at z ^ 2 compared to z 3. 

Finally, we note that the best-fltting scale lengths for 
all our putative Lya halos are relatively small, between 
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TABLE 3 

Large scale Flat-Field Uncertainties 



Field 


Histogram FWHIW 
(lO-^ADU s-i bin-i) 


Median Sky Value 
(ADU s-i pixel-i) 


Uncertainty'' 
(percent) 


Expected Poisson Error 
(10-^ADU s-i bin-i) 


Flux error 
^-|^Q— 19 gj-gg g-i arcsec"^) 


C-03 


1.30 


0.078 


0.167 


0.122 


8.7 


K-03 


1.15 


0.174 


0.066 


0.224 


6.2 


02 


0.75 


0.061 


0.123 


0.087 


9.9 



^ See Figuro[3]and the discussion in the text. 

^ Defined as Histogram FWHM / Median Slcy Value. 
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Fig. 4. — The radial profiles of the LAEs derived from the C-03 sample. For each of these figures, the black circles represent the radial 
profile of the stacked LAEs, the blue open triangles represent the large-scale PSF, as determined in §4, and the red squares show the stacked 
radial profile of the stellar sources found in each frame. LAE bins that have negative flux are denoted by a filled point at the bottom of 
each plot. The large-scale PSF and stellar sources have been normalized to match the three innermost points of the LAE profile. The 
horizontal dashed line represents the flux limit comparable to our systematic uncertainties. The differing profiles relate to the statistically 
complete sub-sample (top left), the UV Bright sub-sample (top right), and the high EW sub-sample (bottom left) and low EW sub-sample 
(bottom right), respectively. 



C/S and I'.'O in angular units, or '-^ 3 to '-^ 8 kpc in phys- 
ical size. This is significantly less than the 20 to 25 kpc 
Lya scale lengths found by Sll in their study of LBGs, 
and only marginally consistent with the 9-20 kpc scale 
lengths found in the LAE samples studied by M12. This 
is easily seen in Figure [7l where the best-fitting exponen- 
tial profiles for our sub-samples are clearly steeper than 
that of Sll and M12. If the ECDF-S LAEs do have ex- 
tended Lya halos, they are smaller than the halos found 



by other studies. 

Could some instrumental or sampling issue be affect- 
ing our conclusions? Since we are directly comparing 
the LAEs to a sample of point sources on the same CCD 
frame, and are using the identical techniques of stacking, 
averaging, and fitting, it seems unlikely that there are ef- 
fects that would cause large systematic errors. However, 
one possible difference between the two samples is in their 
observed flux. The LAEs measured in Figures H] through 
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Fig. 5. — The radial profiles of the LAEs derived from the K-03 sample, similar to Figure|4] 
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TABLE 4 

Point Source Subtracted Exponential Fits 





C-03 


K-03 


02 


Sample 




r,i (kpc) 


a 


r„ (kpc) 


/~i a 


r„ (kpc) 


Statistically Complete 


3.6± 1.1 


5.7 ± 1.0 


3.1 ± 1.5 


4.4 ± 1.0 


1.3 ± 1.6 


5.4 ±3.7 


UV Bright 


3.8 ± 1.0 


6.0 ±0.9 


1.3 ±0.6 


6.4 ± 1.8 


1.5 ± 1.5 


5.7 ±3.2 


High Equivalent Width 


1.5 ± 1.2 


5.5 ± 2.6 


1.3 ±6.8 


2.8 ±5.0 


0.6 ±0.4 


3.7 ± 12 


Low Equivalent Width 


2.4 ± 1.1 


5.8 ± 1.5 


1.2 ±0.3 


8.4 ± 1.5 


0.4 ± 1.6 


5.5 ± 11 



Flu 



; in units of 10 



[S]are, on average, more than 100 times fainter than their 
point source counterparts, and hence have larger posi- 
tional uncertainties. Since the act of stacking relies on 
knowing the objects' centroids, it is possible that in our 
co-addition procedure, we have artificially introduced a 
"halo" into our LAE sample. 

To quantify this potential strength of this effect we per- 
formed a series of Monte-Carlo experiments, each time 
adding 200 faint (S/N = 9) artificial point sources to 
the C-03 frame, and stacking the results, while adding a 
random amount of positional uncertainty, Upixei, into the 
assumed x and y coordinate of each object. We chose 
a range of (Tpixoi values that should span the range of 
any conceivable uncertainty (0.5, 0.75, 1.0, 1.5, and 2.0 



pixels). The fiux of the artificial sources was chosen to 
isolate the effects of centroiding errors from that due to 
signal-to-noise effects. We then performed the identical 
stacking, averaging, sky subtraction, and exponential fit- 
ting steps as before. Any non-zero values of C„ would 
indicate an additional uncertainty above that found from 
the photometric errors. 

As illustrated in Figure [S] and summarized in Table [1] 
this positional broadening has the potential to mimic the 
effects of a diffuse halo: large 1 pixel) uncertainties in 
the {x, y) positions of individual LAEs can cause a point- 
source stacked image to appear extended, and generate a 
false exponential profile. Obviously, the larger the value 
of CTpixci, the larger the size of the apparent halo, so the 
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Fig. 6. — The radial profiles of the LAEs derived from the 02 sample, similar to Figure|4] 
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TABLE 5 

Effects of Centroiding Errors on Exponential Fits to Point 

Sources 



(pixels) 



(10- 



Cn 

'17 —1 

ergs s 



cm 



(arcsec) 



(kpc) 



0.0 

0.5 

0.75 

1.0 

1.5 

2.0 



-0.05 ± 0.32 
0.5 ± 4.7 

1.1 ± 6.7 
1.0 ± 3.6 
0.8 ± 1.7 

1.2 ± 1.8 



1.3 
0.4 
0.4 
0.5 
0.6 
0.6 



6.2 
1.7 
0.9 
0.7 
0.6 
0.4 



11 ± 52 
3.3 ± 14 
3.3 ± 7.3 
4.2 ± 5.8 
5.0 ± 5.1 
5.0 ± 3.6 



The meaning of these parameters are identical to Table [T] 



amplitude of this effect depends strongly on the bright- 
ness distribution of the LAE sample being studied. For- 
tunately, this can be estimated by repeating our artifi- 
cial star experiments us ing LAE magnitude di stributions 
similar to those given in lGronwall et all (|2007f ). and mea- 
suring the uncertainty in their centroid as a function of 
signal-to-noise. The result of this artificial star experi- 
ment is a median coordinate uncertainty of ^ 0.05 pixels 
for objects near the bright end of the LAE luminosity 
function, ^ 0.5 pixels for objects near the faint end of 
the statistically complete LAE sample, and ^ 0.8 pixels 
in the worst case scenario. The curves of Figure |8] then 



suggest that errors in the LAE centroiding can explain 
some, but not all, of the excesses we detect. 

We conclude that we have some, but not conclusive, 
evidence for LAE halos in our sub-samples. LAEs at 
z ^ 3 with low equivalent widths and bright continuum 
magnitudes are more likely to show evidence for an ex- 
tended halo, but due to the uncertainties in flat-fielding, 
sky subtraction, and coordinate centroiding, such results 
are still uncertain. The Lya halos we do detect are much 
smaller in physical scale than those of Sll and M12, up 
to a factor four smaller. There does appear to be clear 
differences between our results at z ~ 3 and z 2, where, 
despite the reduced cosmological effects, we see no evi- 
dence for any extended halo. This may suggest an evo- 
lutionary component to the formation of extended Lya 
halos. 

6. MONTE-CARLO SIMULATIONS 

Figure [3] and Table [T] illustrate the fundamental limita- 
tion of all LAE surface photometry performed over scales 
greater than ~ 100 kpc. It does not address the question 
of small-scale background fluctuations, and it is at least 
conceivable that on galaxy-sized scales, the uniformity 
of the sky background improves. To test this possibility, 
we performed a series of Monte Carlo simulations, using 
model LAE halos of the form suggested by Sll. 

We began by using the GALFIT analysis package 
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Fig. 7. — The residual radial profiles of our stacked LAEs, formed by subtracting a scaled version of each frame's point-source stack from 
the sub-sample's original LAE stack. Sub-samples which display no evidence for extended emission are not shown. Bins with negative fiux 
are denoted by a filled point at the bottom of the plots. The best-fitting exponential model for each sub-sample (given in Table[l{ is shown 
as the solid red line; for comparison purposes, the LAE only exponential fit of Sll (see their Table 2) is overplotted as the green dotted 



line, and the — 1 < (5? 



< 0.5 (see their Table 1) exponential fit of M12 is overplotted as the dashed magenta line. 



(jPeng et al.ll200ai2010f l to create a set of artificial LAEs, 
with each LAE consisting of two components: a point 
source core, and an exponential halo of the form sug- 
gested by Sll (with an adopted scale length, r„ — 
28.4 kpc). In each of the simulation runs, we chose 
a global halo-to-core flux ratio (corresponding to dif- 
fering values of C„), convolved our composite galaxies 
with the frame's large-scale PSF, and assigned these sim- 
ulated LAEs random fluxes drawn from the ISchechteH 
19761) luminosity func tion appropriate for their redshift 
CiarduUo et al.l 720121) . We then randomly placed the 



artificial galaxies on their respective frame and stacked 
those objects brighter than the frame's 90% complete- 
ness limit. We performed eleven such simulations for the 
C-03 and the 02 images, with halo-to-core ratios spaced 
evenly between zero and four. 
The top panels of Figure |9] show the results of these 



simulations for halos of the type suggested by Sll. It 
is clear that at this halo-to-core ratio (corresponding to 
Cn ~ 4.6 X 10^^^ ergs cm~^ s~^ arcsec"^ for the C-03 
stacks, and C„ « 1.1 x 10~^^ ergs cm~^ arcsec"^ for 
the 02 stacks), extended Lya emission would have been 
easily detected in both stacked images and have a mor- 
phology that is quite different from the observed profiles 
displayed in Figures |4][7l The panels also demonstrate 
that the surface brightness limits to which we can confi- 
dently trace Lya emission are similar to those predicted 
by the large-scale flat-fielding error. Even at scales less 
than ^ 100 kpc, the systematic effects discussed in §4 
are still the dominant error terms in our analysis. 

Of course, it is certainly possible that the Lya halos 
of high-redshift LAEs are fainter than those described 
by Sll, and have correspondingly smaller values of C„. 
To examine this possibility, we visually inspected and 
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Fig. 8. — The effects of centroiding errors on our measured pro- 
files. The black filled points represent the co-added profile of 200 
artificial point sources (each with a signal-to-noiso of 9) on the C- 

03 image. The open triangles represent the same data, but with 
differing amounts of Up^^^i , with red, orange, green, blue, and pur- 
ple (running from closest to furthest away from the solid squares) 
corresponding to fJpi^^i = 0.5, 0.75, i.O, 1.5, and 2.0 pixels, respec- 
tively. Each profile has been normalized to the central point for 
comparison purposes. 

repeated our non-linear least squares analysis on the 
simulated halos with smaller halo-to-core ratios. The 
bottom panels of Figure [9] illustrate those models with 
marginally-detected halos. In the C-03 panel, we found a 
limiting value of C„ 2 x 10~^^ ergs cm~^ s~^ arcsec"^, 
while in the 02 plot, the limiting value of C„ is « 

4 X 10~^^ ergs cm~^ s~^ arcsec^^ Any LAE halos with 
fluxes below these limits would be undetectable on our 
images. We note that the stacked Lya halo found by 
Sll for LBGs is relatively bright, as it contains ~ 3.4 
times more flux than the galaxy's central core. The ha- 
los detected by M12, however, are substantially fainter, 
with best fit values for C„ ranging from 3.2 down to 
0.7 X 10~^^ ergs cm~^ s~^ arcsec"^. Therefore, for our 
samples of LAEs, although the bright halos suggested by 
Sll are ruled out by our analysis, fainter halos of the 
form claimed by M12 are still plausible. 

7. DISCUSSION 

Most cosmological simulations - even those which in- 
clude radiative transfer from Lya - do not include the 
effects of resonant scattering of photons from the inter- 
a nd circum-gal actic medium. One exception is the work 
of iZheng et al.l (|201l[ ). who calculated the expected spa- 
tial distribution of extended Lya emission from LAEs at 
z ^ 6. These models predict that in the Lya emission 
line, a high-redshift LAE should have a central cusp, and 
a slowly declining envelope which extends outward from 
^ 0.2 Mpc to '--^ 1 Mpc (in co-moving distance). Unfortu- 
nately, it is difficult to directly compare our observational 
results at 2 < z < 3 to the predictions of these models. 
The distribution of gas at z ~ 6 is largely defined by in- 
fall, and is related to the virial radius of the galaxy; both 
these parameters evolve rapidly in th e Gyr betwe en z ~ 6 
and z ~ 3. In addition, the .Zheng et al.l (|2011|) models 



have a poorer spatial resolution than our data, making it 
difficult to compare except at the largest radius. These 
simulations also do not include the effects of dust, and, 
depending on how this component is distributed, its ef- 
fect on the propagation of L ya photons can be dramatic 
(e.g., iNeufcl d 1991; Finkclst ein et al.ll2009f l. Finally, the 
IZheng et al.l ()2011f ) models do not take into account the 
internal structure and kinematics of the ISM. Thus, it is 
not possible to directly compare their model predictions 
with observational res ults. 

More recently, Diik stra fc Krameil (|2012| ) performed 
radiative transfer calculations for Lya photons prop- 
agating through clumpy, dusty, large scale outflows. 
While their models do predict the existence of ex- 
tended Lya halos, their initial models (I-IH) have ha- 
los that are more concentrated (approximate inner scale 
lengths of ~ 4 kpc) than those observed by Sll, and 
fainter by almost an order of magnitude (« 4.5 x 
lO^'^'^ ergs cm~^ s~^ arcsec^^ at a radius of w 5"). 
More complex models (IV-V) that include the decel- 
eration of clumps are better matches to the halos ob- 
served by Sll, but underpredict the flux by a factor 
of 2-3. Unfortunately, due to our observational limits, 
we would not be able to detect halos of these types; 
only the two innermost points of their halos would 
be measurable in our data, with fluxes of 0.1 — 1 x 
lO"^'^ ergs cm~^ s~^ arcsec"^. There are several the- 
oretical studies of Lya emitters within a cosmological 
framework that also include explicit radiative transfer 
(iForero-Romero et al.|[20TlHGarel et al.ll2012t[Orsi et al.l 
|2012| ) which could be applied to our data, and these 
should assist in the interpretation of our results. 

As described above in detail, accurate surface photom- 
etry at low flux levels is difficult, and fraught with sys- 
tematic errors. Consequently, direct comparisons with 
other Lya halo measurements is difficult and may ig- 
nore imp ortant systematic e ffects. Using long-slit spec- 
troscopy, iRauch et al.l ()2008[ ) found extended Lya emis- 
sion in roughly half of a sample of 27 LAEs in the red- 
shift range 2.7 < z < 3.8. But this emission was weak: 
no brighter than ~ 1.8 x 10~^^ erg s~^ cm~^ arcsec"^ 
at ~ 4" (~ 30 kpc) and well below the level where sys- 
tematic uncertainties dominate our error budget. The 
halos found by M12 using the stacked images of several 
hundred field LAEs identified on Subaru frames are sim- 
ilarly faint: for c ompa rison, the LEG halos of Sll and 
iHavashino et al.l H2004D . which were measured in highly 
overdense regions of the universe, are almost an order 
of magnitude brighter. Other attempts to image LAE 
halos, both from the gr ound ([Smith et al. 2012) , and 
from space (|Bond et al.l 120101: iFinkelstein et al.l i2011| ) 
have placed only weak upper limits on the phenomenon, 
and generally probe very different spatial scales. 

Sll hypothesize that extended Lya halos are a "generic 
property of high redshift star-forming galaxies" , while 
M12 argue that the emission properties of Lya halos de- 
pend on environment. Neither of these studies took into 
account the systematics associated with the large radius 
point-spread-function or large-scale flat fielding, but nei- 
ther are in direct conflict with our results. The median 
LAE in our analysis is over a magnitude fainter than the 
LBGs studied by Sll, and our surface brightness lim- 
its, 9.9 X 10^^^ erg s^^ cm~^ arcsec"^ for z ^ 2.1 and 
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Fig. 9. — A comparison of Monte-Carlo models with an exponential halo added to the large-scale PSF. The meaning of the symbols is 
identical to Figurel4l The top row indicates the halo-to-core ratio consistent with the results of Sll (rn = 28.4 kpc). An extended halo 
of this type would have been clearly detected in our data. The bottom row indicates the models where an extended halo is just barely 
detected, with a halo-to-core ratio of ^ 1.45 in the C-03 plot, and ^ 1.09 in the 02 plot. 



8.7 X 10^^^ and 6.2 x 10^^^ erg s^^ cm^^ arcsec^^ for the 
two z 3.1 samples are brighter than the values quoted 
by M12. Thus, there are several ways to reconcile the 
three surveys. 

First, LAEs may have significantly fainter Lya halos 
than Lyman-break galaxies. In general, the star forma- 
tion rates of LAEs are a n order of magnitude smaller 
than those of LBGs (e.g., IGronwall et al.|[2"007D . and it 
is reasonable to hypothesize that the lower masses as- 
sociated with LAEs also imply less circum-galactic ma- 
terial. These two factors may combine to produce Lya 
halos whose surface brightness is highly sensitive to the 
luminosity of the underlying galaxy. This argument 
is somewhat mitigated by the fact that Lyman-break 
galaxies have more internal extinction than their LAE 
counterparts, but there are ways to distribute this dust 
so that the resonantly scattered Lya photons c an es- 
cape their surroundings (Ncufcld 1991; Finkelstci n et al.l 
l2009f ). The fact that our best evidence for extended emis- 
sion occurs on the C-03 and K-03 stacks of LBG-like ob- 
jects and low-equivalent width sources argues in support 
of this interpretation. 

Alternatively, the observability of Lya halos may be 
due largely to galactic geometry. The probability of de- 



tecting Lya emission from a galaxy depends on many 
factors, including viewing angle, dust content, gas dy- 
namics, and intrinsic geometry. By definition, LAEs are 
systems where the Lya escape path falls along our line- 
of-sight; in galaxies selected via the Lyman-break tech- 
nique, Lya may leak out in other directions, thus fa- 
cilitating the detection of a "halo." Sll's data argue 
against this interpretation, since the halos surrounding 
LBGs with Lya emission are no different from those 
where Lya was only seen in absorption. However, the 
Lya equivalent widths of typical LAEs are much greater 
than those associated with LBGs, so this selection effect 
may still play a role. Recently, Verhamme et al.l ()2012[ ) 
have claimed that Lya emission from galaxies can have 
strong inclination effects, with LAEs with high equiv- 
alent widths being more likely to be found in face-on 
orientations than edge-on orientations. Although there 
are many details still left to be resolved, it is also worth 
noting that the simulations of I Verhamme et al.l (|2012[ ) 
show a substantially weaker halo than Sll. 

A third possibility concerns the environment of the 
LAEs studied. Ml 2 have argued that galactic envi- 
ronment plays an important role in the observability of 
Lya emission. The LBGs studied by Sll are primarily 
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located in overdense region s of the universe, including 
the well-known field SSA22 (jSteidel et al.lll998D and two 
regions centered on hyper-luminous QSOs (HS1700-I-64 
and HS1549-f 195). In these "proto-cluster" environ- 
ments, there is likely to be a large amount of circum- 
galactic material, which could scatter the Lya photons 
which escape from galaxies. Indeed, not only is SSA22 
known for possessing large numbers of luminous Lya 
blobs fe.g.. lMatsuda et al.ll200HSteidel et al.ll2000[ ). but, 
according to M12, its LAEs are similar in strength and 
profile to those found by Sll for the region's LBGs. In 
lower density regions. Ml 2 found the Lya halos of LAEs 
to have steeper profiles with less total luminosity. The 
LAEs studied in this paper are all in the blank-field (com- 
parable to the —1 < (52^^ < 0.5 sub-sample of M12) en- 
vironment of the ECDF-S. If this is true, then the lack 
of halo emission around our galaxies should not be too 
surprising. 

Finally, we note that our Lya halo measurements at 
z 3.1 differ from those at z ~ 2.1: despite the 
(1 + z)** gain in surface brightness, we see no evidence 
for extended Lya halos in any of our lower-redshift sub- 
samples. One possible explanation for this difference 
lies in the uniformity of LAE populations. A number 
of analyses have shown that LAEs si z ^ 2.1 are, on av- 
erage, slightly larger, dustier, and more heterogeneous 
than their z ~ 3.1 counterparts (iNilsson et al.l 120091 
[20Tlt IGuaita et all [MH : IBond et al.l I2012f ). This can 
be understood if the systems that we call LAEs actu- 
ally consists of two separate populations: one containing 
galaxies with low stellar mass and little dust, and the sec- 
ond consisting of objects that are redder, dustier, more 
massive, an d possibly more evoly e d. U nder the inter- 
pretation of IShimizu fc Umemural (|20100, the first type 
of LAEs are young objects undergoing an early burst 
of star formation, while Type II LAEs are produced by 
delayed accretion of sub-galactic halos onto established, 
more massive systems. Consequently, the star formation 
which generates Lya photons in Type II systems occurs 
in the outskirts of systems with lower density, more dif- 
fuse neutral hydrogen halos. If this is the case, then the 
halos of these objec ts would be more difficult to detect. 
Since the models of IShimizu fc Umemural ()2010[ ) predict 
that the fraction of Type II LAE should increase over 
time, their scen ario can explain b oth the evolution of 
LAE properties (|Guaita et al.lf20Tll ) and our lack of halo 
detections. 

8. CONCLUSION 

Surface photometry at very low fiux levels is treach- 
erous, and extreme care must be used in modeling sys- 
tematics such as the large-radius point-spread-function 
and the large-scale flat-fielding errors. By taking such 
factors into account, we were able to form stacked pro- 
files of three sets of LAEs in the Lya emission-line and 
to obtain robust estimates of their uncertainties. At 
z 3.1, we find mixed evidence for the presence of 
extended Lya halos down to a surface brightness limit 
of « 7 X 10~^^ ergs cm~^ s~^ arcsec"^ At z 2.1, 
where cosmological dimming is a factor of 3.2 less, we 
see no extended emission down to surface brightnesses 
of « 1 X 10~^^ ergs cm^^ s^^ arcsec"^. This is the 
first time the rigorous techniques of low-flux level surface 
photometry have been applied to the problem of galactic 



emission-line halos. 

While our observat i ons ar e not in direct conflict wit h 
those of iSteidel et"all (|2011[ ) and lMatsuda et al.l (|2012D . 
they do place important constraints on the systematic 
behavior of Lya halos and imply that such studies have 
higher uncertainties than originally thought. 

However, further progress on this problem will be dif- 
ficult in the near-term. To image the emission-line ha- 
los around LAE and LBG halos, one should obtain deep 
narrow-band exposures on large telescopes with wide- 
field imagers that also have low scattered light profiles 
and very stable flat-fields. This combination is very dif- 
flcult to obtain. 

Alternatively, it is possible to study the emission-line 
halos of high-redshift galaxies with spectroscopy. Start- 
ing in 2014, the integral-fleld units of the Hobby-Eberly 
Telescope Dark Energy Experiment (HETDEX) will ob- 
tain two-dimensional spectra for 1 0^ Lya emitter s 
in the redshift range 1.9 < z < 3.5 (|HiU et al.l[2007[ ). 
These data, which will have ~ 0'.'5 spatial resolution 
and ~ 6 A spectral resolution, should allow us to trace 
circum-galactic Lya emission to levels far below what 
is currently possible and to determine the properties 
of these systems as a function of luminosity, equivalent 
width, dust content, redshift, and environment, and thus 
illuminate material in the high-redshift universe's cosmic 
web. 
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